Scaling properties of pyrex and silicon surfaces blasted with sharp particles 
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The blasting of brittle materials with sharp particles is an important fabrication technology in 
many industrial processes. In particular, for micro-systems, it allows the production of devices with 
feature sizes down to few tens of microns. An important parameter of this process is the surface 
roughness of post-blasted surfaces. In this work the scaling properties of Pyrex glass and silicon 
surfaces after bombardment with alumina particles is investigated. The targets were bombarded at 
normal incidence using alumina particles with two different average sizes, 29/im and 9/im, respec- 
tively. This investigation indicates that the resulting surfaces have multifractal properties. Applying 
multifractal detrended fluctuation analysis (MFDFA) allowed us to determine the singularity spec- 
trum of the surfaces. This spectrum did not depend on the target material or on the size of the 
particles. Several parameters quantifying relevant quantities were determined. We argue that for 
scales below 5/im, fracture processes are dominant while at large scales long range correlations are 
responsible for the multifractal behaviour. 

PACS numbers: 



I. INTRODUCTION 

Powder blasting technology is among several tech- 
niques used in the micro-machining [3, 0, El of devices 
on silicon and other materials. With such abrasive tech- 
niques, one can achieve high erosion rates; higher than 
those that can be obtained with conventional dry, or 
wet etching processes; such as plasma etching or chem- 
ical etching. In the field of micro-electromechanical- 
systems (MEMS), powder blasting has already been 
used for the fabrication of inertial sensors 0, peristaltic 
micro-pumps[3| and miniaturized capillary electrophore- 
sis chips Erosion with sharp particles is also a widely 
used technique in aerospace and automotive industries. 
Because of its involvement in many applications, it is im- 
portant to investigate the surface morphology of blasted 
surfaces resulting from erosion since the performance of 
many devices will depend on the surface roughness. Such 
investigation will also give an insight into the physical 
mechanisms at work during the bombardment of materi- 
als with particles. The mechanisms involved in the ero- 
sion of brittle materials with sharp particles have been 
the subject of several studies. Several models have been 
developed d, 0] that are based on simple static inden- 
tation theory. These models empirically relate the ero- 
sion rate to the material's properties such as the frac- 
ture toughness, the hardness and the Young modulus 
of the material. These models simply state that an in- 
dentation force is generated by the impact of bombard- 
ing particles which, in turn, results in the formation of 
crack patterns. Some cracks penetrate the material ra- 
dially away from the surface into the bulk material; oth- 



ers will nucleate and form a lateral ring parallel to the 
surface 0, H, HI- T he radial cracking process was exten- 
sively studied jiol . [rfl ] and was used as a method for ma- 
terials toughness measurement [12| . Lateral cracks are 
responsible for the removal of material in abrasive and 
wear experiments on brittle materials (l3l [Til [l5| . 
In separate studies, substantial work was dedicated to 
the understanding of fracture surfaces [la ] in brittle ma- 
terials. These studies focus on the scaling properties of 
fracture surfaces resulting from an applied load on the 
material. These surfaces are found to be self-affine, i.e. 
the root-mean-square surface fluctuations, averaged over 
a distance L follows the scaling relation \u\: 



(1) 
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where the scaling exponent a is often called the rough- 
ness or the Hurst exponent. In a variety of materials, 
the value of a was found to be approximately 0.8 over 
two or three decades of scaling range. For this reason, 
it was conjectured to be universal i.e. independent of 
the material (ductile or brittle), the fracture mode and 
the fracture toughness [lj| [HI Il9| . However, this uni- 
versality was questioned since the discovery of a second 
exponent at the nanometer scale [13, HH H^]. The value 
of this second exponent is significantly smaller than 0.8 
and close to 0.5. To explain this fact, it was proposed 
that the fracture front could be imagined as a line mov- 
ing through a random medium. Thus, the evolution of 
the crack front can be described by a local nonlinear 
Langevin equation [2^, [24|, [2|| which predicts a crossover 
between two regimes corresponding to a = 0.5 at small 
scales and to a = 0.75 at large scale. In the framework 
of this model, the crossover length decreases rapidly with 
the crack s pee d which was also predicted by numerical 
simulations [26| . 

The present work is dedicated to the study of the seal- 
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FIG. 1: Surface profiles after blasting of Pyrex and silicon 
with particles of size S = 9 fim (Py9 and Si9), and with parti- 
cles of size 6 = 29 fim (Py29 and Si29). 



ing properties of surfaces resulting from a bombardment 
by sharp particles. Pyrex glass (borosilicate glass) and 
silicon were the materials used in this investigation. We 
will investigate the effect of the material and the size of 
the bombarding particles on those properties. 



Particles size 


9 fim 






29fim 


Pyrex 


0.7 ±0.02 


fim 


2.2 ± 


0.02 fim 


Silicon 


0.65 ±0.01 


/im 


2.3 ± 


0.01 fim 



TABLE I: Values of the rms roughness as a function of the 
particles size and the target material. 



Sample 


a\ 


Q-2 




Si9 


0.77 ± 0.02 


0.4 ± 0.01 


4.8 


Si29 


0.81 ± 0.01 


0.55 ± 0.002 


4.7 


Py9 


0.74 ± 0.02 


0.32 ± 0.003 


4.7 


Py29 


0.8 ± 0.02 


0.53 ± 0.04 


5.2 



TABLE II: Values of the roughness exponents and the 
crossover length £ for the four samples. 



III. SCALING PROPERTIES OF BLASTED 
SURFACES 

A. Global scaling exponents 



II. EXPERIMENTS 

In the experiments carried out, surfaces are exposed to 
a directed particle jet, which results in mechanical ma- 
terial removal. The particles are accelerated towards the 
target with a high-pressure air flow through a circular 
nozzle (with a diameter of 1.5 fim). The particles hit the 
target under normal incidence, with an average speed of 
290 m/s, in a ventilated box. A lateral movement of the 
target ensures an evenly etched surface. The average 
diameter of the bombarding alumina particles was 9/xto 
and 29fim, respectively. We performed measurements 
on one-dimensional cuts of blasted surface using a me- 
chanical surface profiler(Sloan Dektak II), over a length 
of 1 mm. Each scanned profile is made of 8000 data 
points. Typical profiles obtained after blasting Pyrex 
and silicon surfaces by alumina particles are shown in 
figure [TJ Note the difference in the roughness amplitude 
in the two cases corresponding to different particle 
sizes. Large amplitudes are obtained when the targets 
are blasted by large particles. The typical rms surface 
roughness is a = 0.7 fim and a = 0.65/im for Pyrex and 
silicon respectively, when 9fim particles are used. In 
table HI measured values of the surface rms roughness 
are shown, as a function of the bombarding particle 
sizes and the target material. Samples bombarded with 
particles having a diameter S are denoted X 5, where X 
denotes the material's symbol. For example Sz29 means 
silicon bombarded with particles having a diameter of 
29/im. 



We first investigate the global scaling behaviour of the 
blasted surfaces by calculating the height-height correla- 
tion function defined as : 

C(r) = ((h(x + r)-h(x)) 2 ) 1/2 (2) 

where the outer brackets mean the average over all posi- 
tions x and h(x) defines the surface profile. This function 
measures how correlated two points are on the profiles at 
a distance r of each other. If the surface is self-affinc 
then C'(r) oc r Q , over the scaling range. Here, a is the 
roughness exponent to be determined. We computed 
C(r) for one dimensional cuts and averaged the result 
over a total of ten scanned profiles. In figure [2] we plot 
C(r) for 2 fim < r < 60 fim. We notice the existence 
of a cross-over, occurring at a characteristic length scale 
£, separating small and large scales with two different 
roughness exponents, i.e. for r < £, C(r) oc r ai and for 
r > £, C(r) oc r° 2 . The linear fit of log(C) vs log(r) 
in the scaling region, determines the value of the rough- 
ness exponent a\ for small scales and oti at large scale. 
These values are summarized in table [TTl along with the 
characteristic length scale £. 



B. Multifractal properties and singularity 
spectrum 

The global analysis performed above can only reveal 
the existence of two roughness exponents. In many situa- 
tions however, rough profiles exhibit a range of roughness 
exponents. Such profiles are called multifractal while 
profiles exhibiting a single roughness exponent are called 
monofractal|27f. Our aim is to determine the full range 
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FIG. 2: Plot of the logarithm of the height-height correlation 
function versus the logarithm of the distance, for Pyrex glass 
and silicon blasted with alumina particles with sizes of 29 
and 9 microns. This plot reveals the existence of a cross- 
over length scale £ separating two scaling regions having two 
different global roughness exponents. 

of local roughness exponents for the blasted surfaces us- 
ing multifractal analysis. To do so, we characterized 
the profiles by computing the so-called singularity spec- 
trum, which determines the distribution of the whole 
range of local roughness exponents. We used the mul- 
tifractal detrended fluctuation analysis (MFDFA) [28| . 
This method has become popular thanks to its simplic- 
ity and its easy computer implementation. Other meth- 
ods exist such as the wavelet transform modulus max- 
ima methodfWTMM) [H but the MFDFA has several 
advantages [30] . The MFDFA is the extension of the de- 
trended fluctuation(DFA) [3l[ method which was used 
to compute the roughness exponent of monofractal sig- 
nals and for the identification of long range correlations 
in non-stationary time series [28| . MFDFA is an efficient 
tool to eliminate undesirable trends in fluctuations. This 
method applied to our experimental data can be sum- 
marized as follow: Given a profile h(i), i = 1,...N, we 
compute the integrated profile, 

N 

H( l ) = j2my<h>) (3) 

i=l 

where N is the number of data points and < h > is 
the mean height of the profile. The the whole profile is 
subdivided into Ml — Nj L non-overlapping segments of 
length L (here L is the number of data points in each 
segment). Since Ml is not always an integer, some data 
will be ignored during this procedure. To take them into 
account, the subdivision is performed from both ends of 
the profile, which results into 2Ml segments. In each 
segment n the polynomial trend P n is subtracted from 
the data. This polynomial is determined by a least square 



fit to the data in each segment. Polynomials of degrees 
higher than 1 can be used, corresponding to MFDFA2, 
MFDFA3, etc. After detrending in each segment, the 
variance of the result is calculated : 

1 3=L 

F\n, L) = i^[ff((n - 1)L + j) - P n (j)} 2 (4) 
i=i 

This expression is then averaged over all segments n and 
the value of the q'^-order fluctuation function is calcu- 
lated [||: 

Here q is a real number. For a fractal profile F q (L) follows 
a power law relation at large scales i.e.: 

F q (L) ~ L h ^ (6) 

The exponent h(q) is called the "generalised Hurst ex- 
ponent" . For a monofractal profile, h(q) = const 
while for a multifractal profile h is a function of q. For 
positive values of q, h(q) describes the scaling behaviour 
of the segments with large fluctuations, while for nega- 
tive values of q, h(q) describes the scaling behaviour of 
the segments with small fluctuations. The singularity 
spectrum fja) is calculated by performing the Legendre 
Transform [32j of (q, r(g)) with r(q) — qh(q) — l, resulting 
in: 

dr(g) 
dq 

f(a) = qa-r(q) (7) 

The spectrum f(a) can be interpreted as the fractal di- 
mension of a subset of points in the profile characterized 
by the singularity strength a (the local roughness expo- 
nent). For a monofractal profile, a = ao and /(a) = 1, 
where ao is the roughness exponent of the profile. The 
strength of the multifractality of a profile can be char- 
acterized by the difference between the maximum and 
minimum values of a, i.e. Aa = a max — a m i n . 
Figure [3] shows the logarithm of the fluctuation function 
(F q (d)) versus the logarithm of the distance d — LAx, 
where Ax is the spatial increment of the profile, for Py29 
sample. Also shown is the linear regression fit at large 
scale for each value of q. We used a detrended polyno- 
mial of degree 1( MFDFA1), but the result remains un- 
changed when polynomial of degree two and three were 
used(MFDFA2,MFDFA3). We can see clearly that at 
large scale, the fluctuation functions F q (L) are straight 
lines in the double logarithmic plot, with different slopes 
indicating the presence of multifractality. In figured] the 
plot of r(q) versus q is shown for the four samples Py9, 
Py29, Si9 and Si29. Note that for \q\ > 5, the function 
r(q) coincides with its asymptotic form which is a linear 
function of q (45[. Hence, we choose — 5 < q < 5. We no- 
tice clearly that r(q) is a concave function of function of 
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FIG. 3: The logarithm of the fluctuation function F q (LAx) 
as a function of the logarithm of the distance d — LAx, for 
Py29 sample, for q=-5, -3, -1, 1, 3 and 5. The straight lines 
are the regression fits to the data at large scale. The slopes 
of the lines determine h(q). 
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FIG. 5: The singularity spectrum for the four samples. In 
each case, the spectrum was averaged over 12 profiles con- 
taining 8000 points. 
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FIG. 4: The q dependence of r(q) = qh(q) — 1 for the four 
samples. 



q, typical of multifractal profiles. From equations ((7|) we 
can estimate the singularity spectrum for the four sam- 
ples as shown in figure [3 This spectrum is averaged over 
12 profile containing 8000 points each. In table ITTTl we 
summarize the values of the strength of the multifractal- 
ity Aa for the four samples. 



Sample 


Py9 


Si9 


Py29 


Si29 


Aa 


0.94 


0.8 


0.91 


0.94 



TABLE III: Values of the multifractality strength Aa 
Q-min — Q-max for different samples. 



IV. DISCUSSION 

The values of the two global roughness exponents were 
determined for all four samples Py9, Py29, Si9 and Si29 
using the height-height correlation function. All samples 
show a crossover behavior corresponding to two values 
of the roughness exponent a\ and «2 shown in table [U 
The value ol\ is, within the error range, independent of 
the material or the size of impacting particles. For all an- 
alyzed samples, the small scale value is consistent with 
the approximate universal value of 0.8 found for three 
dimensional fracture of brittle materials fl8l l33l |3~H . This 
provides evidence that at scales below the characteris- 
tic length £ ~ 5^m, the dominant mechanism in powder 
blasting is the fracture formation. This observation is 
in agreement with the static indentation theory 0, i, [| 
where the dynamic load of impacting particles creates a 
local load which increases the local stress resulting in the 
formation of lateral cracks, which are responsible for the 
material removal. In general, the universal value of the 
roughness exponent 0.8 corresponding to three dimen- 
sional fracture surfaces is observed at larg e length scales 
(from 0.1 fim to 1 mm, see reference 1 1 (ij : . In contrast, 
the dynamic load of impacting particles induces a cross 
over to smaller values of the global roughness exponent, 
above 5 /zm( see table lITf which are 0.4, 0.32, 0.55 and 
0.53 for Si9, Py9, Si29 and Py29, respectively. The large 
scale values of the global roughness exponent obtained 
for Py9 and Si9 are smaller than those obtained for Si29 
and Py29. The effect of impacting particle's size is felt at 
large scales. Larger particles give larger roughness expo- 
nents, independent of the material being Silicon or Pyrex 
glass. 

We also showed that the surface generated by particles 
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FIG. 6: Top plot, the surface profile of Py29 sample with 
the corresponding local gradient (bottom plot) showing large 
fluctuations of the latter which constitutes the signature of 
intermittency. 

blasting share a common property of multifractality. The 
MFDFA method was used to uncover this property and to 
determine the singularity spectrum for each four samples. 
The strengths of multifractality represented by Aa, as 
displayed in table lllll are close to each other regardless of 
the size of the particles and the material. This multifrac- 
tality could be interpreted in terms of spatial intermit- 
tency. This concept was argued by Krug [35[ to describe 
the scaling of surfaces generated by epitaxial growth 
models incorporating very limited atomic mobility lead- 
ing to a violent spatial intermittent effects and multifrac- 
tal surfaces. This description was borrowed from fluid 
turbulence owing the similarity between Galilean invari- 
ance of turbulent fluids and translational invariance of 
interfaces (36|. In our case, we highlight large fluctuations 
by considering the step size or the gradient at the posi- 
tion x,g(x) = \dh/dx\ for a profile h{x). This quantity 
is the analog of energy dissipation e = (dv/dx) 2 , where 
v is the local velocity of a turbulent fluid [371] . Figure [5] 
shows the local gradient g(x) for Pyrex glass bombarded 
by 29/im-size particles, showing large gradient fluctua- 
tions which suggests spatial intermittency. In addition, 
we consider the distribution of the local gradient in anal- 
ogy with local velocity gradient in fully developed tur- 
bulence, which is described by a stretched exponential 
distribution [38l. [39|. A very similar behaviour is charac- 
teristic of all four samples Py9, Py29, Si9 and Si29 as 
shown in figure where the local gradient distribution 
fits very well with the stretched exponential function: 

P(<?) = iexp(-a<f) (8) 

The fit of the local gradient distribution to equation © 
gives the values of the stretching exponent 7, which are 
1.41±0.05, 1.23±0.03, 1.37±0.05 and 1.21±0.03 for Si9, 
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FIG. 7: Local gradient distribution of the four samples 
Py9(D), Py29(0), Si9(0) and Si29(A). Continuous lines are 
the fit to the stretched exponential distribution (|8]). 



Py9, Si29 and Py29 respectively. The form of the local 
gradient distribution suggests that non-linearities must 
be present at large scale, leading to h — > —h symmetry 
breaking. Indeed, the skewness s of the profiles for all 
samples are non-zero and have the values s ~ -0.4, -0.35, 
-0.28 and -0.4 for Si9, Si29, Py9 and Py29 respectively. 
We can now ask the question: What is the origin of the 
observed multifractality? It well known that multiplica- 
tive cascades [37], Eo[ models generate processes known 
to have intermittent behaviour and a multifractal char- 
acter, mirroring the presence of intermittent fluctuations 
with long-range correlations. These long-range correla- 
tions are generated at large scale by features that hierar- 
chically cascade their influence to smaller scales. To de- 
tect the presence of long range correlations in our surface 
profiles, we perform the following testjU|: we generated 
a surrogate data set by shuffling the height data in each 
profile. The newly generated data set preserves the distri- 
bution of the height but destroys the long range correla- 
tions, which means that the surrogate profiles will exhibit 
a monofractal behavior, if the multifractality originates 
from the long-range correlations and not from the height 
distribution. We performed the MFDFA and found that 
the surrogate data set is monofractal with the rough- 
ness exponent a = 0.5 for all for samples. Thus, the 
observed multifractal behaviour is a result of long range 
correlations since the shuffling procedure preserves the 
height distribution. The effect of long-range elastic inter- 
action in developing long-range correlations was reported 
by some authors in the case of elastic chains driven in a 
quenched random pinning [42[, or during crack propa- 
gation [43|| . Undoubtedly a detailed theoretical investi- 
gation is needed in order to determine the origin of the 
long-range correlations in sharp particle's bombardment 
of brittle materials. 
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V. CONCLUSION 

In conclusion we conducted a detailed scaling analy- 
sis of surfaces of two brittle materials, Silicon and Pyrex 
glass, after bombardment with alumina particles of two 
different sizes and 29/xm. The bombardment results 
in multifractal surfaces. This multifractality is common 
to all samples regardless of the nature of the material 
or the size of the particles. We determined the corre- 
sponding singularity spectrum revealing a broad range of 



scaling exponents. We argued that for scales below 5/xm, 
fracture processes are dominant while at large scales long 
range correlations are responsible for the observed mul- 
tifractal behavior. 
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